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This  paper  presents  a  detailed  algorithm  to  compute  a  mean  visibility  report  between  a  range-restricted  ground 
station  and  a  space  object,  generated  by  the  method  of  ratios.  This  method  is  based  upon  the  limiting  orbital 
characteristics  of  a  space  object  versus  a  ground  station,  which  are  used  to  generate  a  visibility  report  from  an 
analytical  model.  For  this  study,  the  visibility  report  is  defined  as  the  mean  minutes-per-day  the  space  object  is 
above  the  site’s  horizon  and  within  sensor  range.  The  computer  simulation  models  orbital  motion  using  first 
order  secular  perturbations  caused  by  mass  anomalies.  The  method  of  ratios  reduces  computation  time  by  over 
97%  when  the  results  are  compared  to  a  visibility  truth  table  created  using  the  satellite  rise-set  algorithm 
developed  by  Alfano  et  al.  This  method  can  also  be  used  to  evaluate  all  objects  in  a  space  catalogue  against  a 
network  of  range-restricted  ground  stations.  Operations  personnel  can  then  use  this  master  visibility  report  to 
maintain  a  current  database  of  all  space  objects  visible  to  a  particular  ground  station,  which  is  central  for  space 
surveillance  radar  sensors  to  perform  efficient  tracking. 
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The  Method  of  Ratios 

David  Negron,  Jr.,1  Salvatore  Alfano,2  and  Daniel  Donald  Wright  III3 


Abstract 


This  paper  presents  a  detailed  algorithm  to  compute  a  mean  visibility  report  between  a 
range-restricted  ground  station  and  a  space  object,  generated  by  the  method  of  ratios. 
This  method  is  based  upon  the  limiting  orbital  characteristics  of  a  space  object  versus  a 
ground  station,  which  arc  used  to  generate  a  visibility  report  from  an  analytical  model. 
For  this  study,  the  visibility  report  is  defined  as  the  mean  minutes-per-day  the  space  ob¬ 
ject  is  above  the  site's  horizon  and  within  sensor  range.  The  computer  simulation  models 
orbital  motion  using  first  order  secular  perturbations  caused  by  mass  anomalies.  The 
method  of  ratios  reduces  computation  time  by  over  97 9c  when  the  results  are  compared  to 
a  visibility  irulh  table  created  using  the  satellite  rise-set  algorithm  developed  by  Alfano 
et  al.  This  method  can  also  be  used  to  evaluate  all  objects  in  a  space  catalogue  against  a 
network  of  range-restricted  ground  stations.  Operations  personnel  can  then  use  this  mas¬ 
ter  visibility  report  to  maintain  a  current  database  of  all  space  objects  visible  to  a  particu¬ 
lar  ground  station,  which  is  central  for  space  surveillance  radar  sensors  to  perform 
efficient  tracking. 

Introduction 


Space  operations  personnel  performing  space  track  have  enjoyed  the  luxury  of 
working  with  a  small  database  of  only  a  few  thousand  space  objects.  However, 
international  interest  in  space  and  the  increased  domestic  reliance  on  orbital  re¬ 
mote  sensing  and  communications  has  caused  the  population  of  artificial  Earth 
satellites,  which  is  continuously  growing,  to  total  more  than  seven  thousand  ob¬ 
jects.  This  growth  of  orbiting  objects  requires  a  revised  strategy  in  tasking  range- 
restricted  radar  sensors  and  maintaining  their  orbital  element  database.  This 
becomes  tremendously  important  for  sensor  sites  with  the  dual  responsibility  of 
space  track  and  missile  warning.  Each  time  a  site  with  these  responsibilities 
identifies  a  space  object  with  elements  that  do  not  match  a  member  of  its  data- 
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base,  operations  personnel  must  transmit  a  message  to  alert  Cheyenne  Mountain 
AFB  orbital  analysts  of  the  uncorrelated  contact.  If  the  site’s  database  is  not 
current  this  becomes  a  frequent  task. 

Sensor  tasking  and  database  maintenance  require  the  processing  of  each  object 
in  a  space  catalogue  against  a  network  of  sites.  The  processing  algorithm  must  be 
both  efficient  and  robust  because  a  catalogue  can  include  thousands  of  objects 
spanning  the  range  of  closed  orbits,  as  shown  in  Figs.  1  through  3.  One  possible 
method  to  create  a  visibility  report  is  to  use  a  step-by-step  propagation  scheme  to 
simulate  a  sufficiently  large  number  of  passes  of  a  space  object  versus  a  station 
and  record  statistics  about  the  angular  separation  of  their  Earth  centered  inertial 
position  vectors.  A  similar  method  involves  creating  a  report  of  topocentric  slant 
range  and  azimuth-elevation  angles  to  determine  if  and  when  an  object  becomes 
visible  to  a  site.  Two  drawbacks  common  to  these  methods  are  lengthy  computa¬ 
tion  time  and  the  dependence  on  periapsis  rate  to  arrive  at  a  reasonable  answer. 
As  periapsis  rate  tends  to  zero,  the  simulation  runtime  becomes  prohibitively 
long  because  the  accuracy  of  the  visibility  report  depends  on  modeling  all  possi¬ 
ble  passes  of  an  object  versus  a  stat  on.  A  geometric  solution  to  determine  if 
there  is  visibility  between  a  range-restricted  station  and  a  space  object  involves 
creating  a  right  triangle  at  the  site  with  hypotenuse  equal  to  the  objects  maxi¬ 
mum  radial  distance  from  the  geocenter.  If  the  vertical  member  of  the  triangle 
becomes  greater  than  one  Earth  radii,  the  space  object  will  have  a  positive 
topocentric  elevation.  This  method  gives  no  insight  about  contact  duration, 
especially  for  a  range-restricted  site. 


Inclination  (degrees) 

FIG.  2.  Space  Object  Distribution  versus  Orbital  Inclination. 


Mean  Motion  (rev/solar  day) 

FIG.  3.  Space  Object  Distribution  versus  Orbital  Mean  Motion 
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The  method  used  to  determine  sensor  tasking  and  database  maintenance  must 
minimize  simulation  runtime  while  maximizing  the  accuracy  of  the  site-object 
visibility  report.  The  number  of  artificial  Earth  satellites,  which  is  continuously 
growing,  renders  the  multiple  pass  and  topocentric  range-angle  simulation  an  ex¬ 
pensive  choice.  The  published  literature  about  space  sensor  and  land  mass  line- 
of-sight  has  centered  about  satellite  deployment  strategies  to  optimize  remote 
sensing  parameters  for  local  and  global  Earth  coverage.  Examples  of  this  work 
include  [1-4]  which  address  how  to  create  optimal  geometry  for  visibility  be¬ 
tween  a  spaceborne  sensor  and  an  Earth  target.  This  paper  addresses  the  oppo¬ 
site  problem  because  the  sensing  platform  is  now  Earth  located  and  its  sensor 
range  is  limited.  The  objective  of  this  paper  is  to  present  an  analytical  model 
that  accomplishes  sensor  tasking  and  database  maintenance  efficiently  enough 
to  execute  on  a  personal  computer. 

The  Ratio  Equation 

Presented  below  is  a  geometric  model  to  formulate  the  ratio  equation. 
Hayes  [1,2]  used  the  concept  of  ratios  to  analytically  study  the  contact  between 
a  spaceborne  sensor  and  an  Earth  target.  What  follows  is  an  independent  deriva¬ 
tion  of  this  equation  and  a  detailed  description  of  its  use  to  detemine  the  visibil¬ 
ity  between  a  range-restricted  site  and  a  space  object.  The  fundamental  unit  of 
length  for  the  ensuing  equations  is  one  Earth  radii. 

The  method  of  ratios  is  based  on  the  assumption  that  the  space  object  is  visi¬ 
ble  to  all  sites  with  latitude  L  and  sensor  range  pnA.x  with  similar  frequency.  Con¬ 
sidering  the  Earth  as  a  fixed  sphere,  the  space  object  is  advanced  in  its  orbital 
plane  with  visibility  from  the  range-restricted  site  being  evaluated  at  n  discrete 
points.  The  mean  anomaly  is  chosen  as  the  controlling  variable  for  this  algorithm 
because  a  constant  step  size  allows  the  object  to  linger  over  a  hemisphere  as  de¬ 
termined  by  its  orbital  eccentricity  and  periapsis.  To  account  for  apsidal  drift, 
the  argument  of  periapsis  is  initialized  in  the  southern  hemisphere  and  the  visi¬ 
bility  contact  determined  for  one  orbit.  This  process  is  repeated,  moving  periap¬ 
sis  in  m  discrete  increments  until  completing  one  cycle.  The  mean  anomaly  and 
periapsis  step  sizes  are  designer  chosen;  for  this  study  n  is  1(30  and  m  is  8.  Ad¬ 
vancing  periapsis  by  45°  results  in  symmetry  between  the  range-restricted  site 
and  the  space  object,  with  only  five  iterations  needed  to  complete  the  simulation. 

As  shown  in  Fig.  4,  the  (fi.fr.fj)  coordinate  system  is  used  to  express  the  ge¬ 
ometry  between  a  site  with  sensor  range  pmx  and  a  space  object.  This  coor¬ 
dinate  system  is  created  such  that  its  principal  axis  is  aligned  on  the  sensor 
meridian,  the  equator  serves  as  the  fundamental  plane,  and  the  origin  is  at  the 
geocenter.  Thus,  from  Fig.  4,  the  sensor  site  and  space  object  position  vectors 
are  given  by 

R  sin.  =  cos(/L)fi  +  sin  (L){s,  (1) 


and 


R  =  R  cos(L')cos  -  (\  +  R  cos(T’)sin  -]f2  +  R  sin(L'){u 
V  2/  \2I 


(2) 
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FIG.  4.  Geometry  Between  a  Range-Restricted  Sensor  Site  and  a  Space  Object  Expressed  in 
the  (fi,  (:.(})  Coordinate  System. 


where  L  and  L'  are  the  latitudes  ot  the  site  and  object,  respectively,  and  8/2  is 
the  angular  measure  between  R  and  R.«r*  in  the  fundamental  plane.  From  the 
inner  product,  the  angle  A  is  given  by 


cos(A)  = 


cos(L)cos(L')cos 


+  sin(L)sin(L'). 


(3) 


Also  from  the  law  of  cosines 


cos(A)  = 


+  1  ~  PM.-t.V~ 

2  R 


(4) 


Equating  these  two  expressions  for  cos(A)  and  solving  for  8  yields  the  ratio 
equation 

e  =  2cos  ’(^-  +  1  -  -  sin(L)sin(L')\ 

\  cos(L)cos(Z.')  / 

The  above  equation  is  valid  when  the  object  is  above  the  site’s  horizon  and 
within  sensor  range,  as  defined  by  the  constraints 

Rcos(L  -  L')  >  1.0,  (6) 


and 


R2  +  1  ~  2Rcos(L  -  I.  )  £  pmx2- 


(7) 


302 


Negron,  Alfa  no,  and  Wright 


For  a  site  with  a  sufficiently  large  sensor  range  viewing  a  low-Earth  orbiting 
space  object,  it  is  possible  for  the  object  to  go  below  the  horizon  prior  to  exceed¬ 
ing  the  site's  maximum  sensor  range.  Should  this  occur,  the  constraint 

P max'  >  (R2  -  1)  (8) 


is  satisfied  and  the  ratio  equation  becomes 

_J\/R  ~  sin(L)sin(L") 
\  cos(L)cos(E’) 


8  =  2  cos 


(9) 


Since  the  object  is  assumed  to  be  visible  with  similar  frequency,  it's  likelihood  of 
being  seen  by  a  specified  site  is  given  by 


/*(  visibility)  = 


2  -rr 


(10) 


The  ratio  equation  presents  numerical  difficulties  when  the  absolute  value  of 
the  inverse  cosine  argument  becomes  greater  than  one,  which  can  occur  fre¬ 
quently  for  a  site  with  infinite  sensor  range  placed  near  a  pole.  To  avoid  this  con¬ 
dition,  it  is  necessary  to  understand  how  the  argument  represents  the  physical 
relationship  between  the  range-restricted  site  and  the  space  object.  If  the  in¬ 
verse  cosine  argument  becomes  greater  than  or  equal  to  one,  then  the  object  is 
completely  out  of  view  and  0  is  zero.  If  the  argument  becomes  less  than  or  equal 
to  negative  one,  then  any  site  at  that  latitude  can  view  the  object  and  0  be¬ 
comes  2v. 


Visibility  Filter 

An  algorithm  to  filter  space  objects  that  are  never  visible  to  a  sensor  site  with 
latitude  L  and  range  pMAX  is  needed  to  eliminate  unnecessary  computer  process¬ 
ing.  For  a  visibility  report  to  be  nonzero,  the  space  object  must  (1)  have  a  posi¬ 
tive  topocentric  elevation  at  some  point  in  its  orbit,  and  (2)  be  within  sensor 
range.  Equations  for  prescreening  each  of  these  constraints  are  presented  below. 

To  determine  if  an  object  has  a  positive  topocentric  elevation,  compute  the 
minimum  geocentric  range, 

Rmis  =  sec  (a),  (II) 


where 


and 


\L\  —  L  MAX  \L\  >  L  MAX 
0  otherwise 


f  inclination 


direct  orbit 


(12) 


(13) 


MAX  |  ir  —  inclination  indirect  orbit 
Also,  the  maximum  geocentric  range,  RM<x ,  for  an  orbit  with  L'max  is  given  by 

Rmax  — 


cos(a)  +  Vp„„5  -  sin5(a)  pMAx  ^  sin2(a) 
0  otherwise 


(14) 
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Thus,  if 


a(l  +  e)  <  R mik 


(15) 


or 


a(l  -  e)  >  Rmax,  (16) 

where  a  is  the  orbital  semi-major  axis  and  e  is  the  eccentricity,  then  the  object  is 
never  visible  to  the  site,  resulting  in  a  visibility  report  of  zero  contact  time. 


Visibility  Criterion 


As  stated  previously,  the  method  of  ratios  is  based  on  the  assumption  that  the 
space  object  is  visible  to  all  sites  with  latitude  L  and  sensor  range  pmax  with  simi¬ 
lar  frequency.  To  determine  if  this  assumption  is  met,  a  visibility  criterion  is  de¬ 
veloped  below. 

An  object  is  visible  to  all  sites  on  a  specific  latitude  with  similar  frequency  if  it 
is  seen  by  the  ith  and  yth  sites  before  the  groundtrace  pattern  repeats.  Define 
0MAX  as  the  maximum  value  of  all  discrete  values  of  8  which  are  available  from 
computing  the  ratio  equation.  If  O.max  is  2tr,  then  at  some  point  in  its  orbit  the 
object  is  visible  to  all  sites  with  latitude  L  and  sensor  range  p* IAX;  hence  the  vis¬ 
ibility  criterion  is  satisfied.  For  8 x  less  than  2tt  the  minimum  number  of  orbital 
revolutions  for  complete  latitudinal  viewing  becomes 


2v 

&MAX 


(17) 


To  determine  if  the  groundtrace  pattern  repeats  before  obtaining  complete 
latitudinal  viewing,  it  is  necessary  to  compute  the  absolute  movement  of  the 
ascending  node  per  orbital  revolution,  which  is  given  by 

n*  =  |»«  -  n|/P.  (is) 

where  tvq,  is  the  rotation  rate  of  the  Earth,  17  is  the  ascending  node  rate  in 
inertial  space  caused  by  the  Earth’s  oblateness,  and  IP  is  the  orbital  period. 
Rounding  to  the  nearest  integer  gives  the  minimum  number  of  revolu¬ 

tions  of  the  ascending  node  about  a  rotating  earth,  Mi*.  to  ensure  complete 
latitudinal  viewing.  Numerically  this  is  determined  from  the  equation 

x,  =  ~,  j  -  1,2 . N,^.  (19) 

u.n 


where  if  any  x,  is  an  integer  in  the  range  0  to  exclusive,  then  the  orbital 

trace  repeats  too  soon.  An  algorithm  for  the  visibility  criterion  is  given  in 
Appendix  B. 


Simulation  Results 


For  this  study,  the  truth  data  is  generated  by  the  method  of  parabolic  blending 
developed  by  Alfano  et  al,  [5],  to  determine  the  rise-set  times  of  satellite-ground 
station  visibility  periods.  The  time  step  used  for  the  blending  is  250  seconds  with 
the  object  and  site  initialized  at  the  vernal  equinox.  The  orbital  simulation  to 
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generate  the  truth  data  includes  first  order  secular  perturbations  caused  by  mass 
anomalies  [6]  where  the  rate  terms  are: 


where  n  is  the  anomalistic  mean  motion,  nu  is  the  mean  motion  at  epoch,  J:  is  the 
second  harmonic coefficient,  p  is  the  semi-iatus  rectum,  e  is  the  eccentricity,  /  is 
the  inclination.  (1  is  the  nodal  rate,  and  w  is  the  periapsis  rate. 

Classical  orbital  elements  from  the  United  States  Space  Command  space  object 
catalogue  are  used  as  test  data  for  this  study,  and  are  summarized  in  Table  1. 
This  subset  of  classical  elements  is  derived  by  sorting  the  catalogue,  which 
includes  seven  thousand  objects,  by  eccentricity,  inclination,  and  mean  motion. 
For  completeness  the  first,  central,  and  last  orbital  element  set  from  each  sort 
result  is  then  processed  against  the  sites  listed  in  Tables  2  through  5.  The  sensor 
site  latitudes  and  range  restrictions  for  the  test  data  are  designer  chosen  to 
demonstrate  the  effectiveness  and  versatility  of  the  method  of  ratios;  hence, 
these  sites  are  fictitious.  This  method  reduced  computation  time  by  97.65% 
when  compared  to  parabolic  blending  for  a  three  month  interval;  a  time  savings 
of  99.41%  was  realized  if  the  truth  data  was  computed  for  an  entire  year. 

To  assist  operations  personnel  with  interpreting  the  method  of  ratios  results,  a 
column  entitled  “Confid.-re”  >«  included  ,vi*h  each  visibility  report.  Each  entry 
in  the  confidence  column  is  made  up  of  a  set  of  flags  indicating  the  results  ol 
three  tests:  (1)  visibility,  (2)  critical  inclination,  and  (3)  non-geostationary  orbit. 
These  tests  serve  to  alert  the  operator  of  repeating,  frozen,  and  geostationary 
orbits,  respectively.  Because  satellites  are  deployed  with  only  finite  precision,  the 
criteria  to  declare  an  orbit  inertially  or  Earth-fixed  are  designer  chosen.  For  this 

TABLE  I.  Classical  Orbital  Elements  (w  =  0  =  MA  =  0°) 


Object 

Mrev/solar  day) 

c 

t(deg) 

I 

1.00272141 

.0000032 

0.0956 

2 

8.36589235 

.0080158 

90.0175 

3 

0.24891961 

.9363060 

64.9874 

4 

0.21467209 

.0668128 

57.35IX) 

5 

13.37659679 

.0145072 

90.2619 

6 

16.09769232 

.0078742 

82.8709 

7 

1.00271920 

.0003109 

0.0099 

8 

12.41552416 

.0036498 

74.0186 

9 

13.84150848 

.0048964 

144.6414 
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TABLE  2.  Visibility  Report  (min/day)  for  Site  at  80°  Latitude  with  1,000  KM  Range 


Space  Object 

0.5  years 

1  (1  year 

5  years 

Ratio  method 

Confidence 

i 

0.0 

0.0 

0.0 

0.0 

— 

2 

01) 

0.0 

0.0 

0.0 

— 

3 

0.0 

0.0 

0.0 

0  0 

— 

4 

0.0 

0.0 

0.0 

0.0 

— 

5 

0.0 

0.0 

0.0 

0.0 

— 

6 

32.1 

32.0 

32.2 

32  d 

Y  Y  Y 

7 

0.0 

0.0 

0.0 

0.0 

— 

H 

0.0 

0.0 

0.0 

0.0 

— 

d 

0.0 

0.0 

0.0 

0.0 

— 

TABLE  3.  Visibility  Report  <min/dav)  for  Site  at  40°  Latitude  with  2,000  KM  Range 


Space  Object 

0  5  years 

1.0  year 

5  years 

Ratio  Method 

Confidence 

i 

0.0 

0.0 

0.0 

0.0 

— 

■> 

0.0 

0.0 

0.0 

0.0 

— 

3 

0.0 

0.0 

mi 

0.0 

— 

4 

0.0 

0.0 

0.0 

0  0 

— 

5 

17.4 

17.4 

17.4 

1 7  4 

YYY 

6 

23.2 

23  1 

23.2 

2.3.1) 

YYY 

7 

0.0 

0.0 

0.0 

0.0 

— 

8 

11  1 

111 

hi 

II. 3 

Y  Y  Y 

d 

27.5 

27.5 

27.5 

27.3 

YYY 

TABLE  4.  Visibility  Report  (min/day)  for  Site  at  0°  Latitude  with  4.000  KM  Range 


Space  Object 

0.5  years 

1 .0  year 

5  years 

Ralio  Method 

Confidence 

1 

0.0 

0.0 

0.(1 

0.0 

— 

2 

3.3 

3  .3 

3.4 

3.5 

YYY 

3 

0.3 

0.3 

0.3 

0.3 

YNY 

4 

0.0 

0.0 

<1.0 

0.0 

— 

5 

hd.H 

60  X 

6d.K 

60.5 

YYY 

6 

17  4 

17.3 

17.2 

17.5 

YYY 

7 

0.0 

0.0 

0.0 

0.0 

— 

8 

68.0 

68.0 

68.0 

67.8 

YYY 

d 

II7..3 

117.4 

117.4 

116.6 

YYY 

study,  an  orbit  with  an  inclination  within  two  degrees  of  the  critical  inclination 
and  an  eccentricity  greater  than  0.1  fails  the  second  test,  an  example  being  space 
object  #3.  Also,  an  orbit  fails  the  non-geostationary  test  if  it  has  an  inclination 
less  than  10  degrees,  an  eccentricity  less  than  0.001,  and  a  period  within 
15  seconds  of  the  sidereal  period.  Processing  the  Space  Command  catalogue 
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TidLE  5.  Visibility  Report  (min/day)  for  Site  at  -60°  Latitude  with  1,000,000  KM  Range 


Space  Object 

0.5  years 

1.0  year 

5  years 

Ratio  Method 

Confidence 

1 

1440.0 

1440.0 

722.9 

574.1 

NYN 

2 

379.5 

584.3 

384  2 

384.5 

Y  Y  Y 

3 

656.0 

062.2 

643.7 

687.5 

YNY 

4 

6K5.9 

685.4 

687.5 

687.1 

Y  Y  Y 

5 

17K. 0 

182.0 

183.2 

183  0 

Y  Y  Y 

6 

.47  6 

37.7 

37.7 

37.4 

Y  Y  Y 

7 

1440.0 

1440.0 

720.8 

574. 1 

NYN 

8 

IKK.K 

188.8 

189.5 

189.6 

YYY 

9 

25.1 

24.0 

24.9 

25.0 

Y  Y  Y 

against  these  two  tests  indicate  that  five  percent  of  the  catalogued  objects  fail 
the  critical  inclination  test  and  two  percent  fail  the  non-geostat ionary  test.  It  is 
worthy  to  note  that  a  circular  or  inertially  fixed  orbit  requires  the  processing  of 
only  one  orbital  revolution  to  produce  an  accurate  visibility  report.  The  reports 
generated  by  the  method  of  ratios  and  the  multiple  pass  simulation  are  nearly 
identical  for  a  space  object  that  passes  all  three  tests,  indicated  by  YYY  in 
Tables  2  through  5.  Should  the  visibility  prescreen  determine  zero  contact,  a  field 
of  dashes  appears  in  the  confidence  column. 


Closing  Remarks 

This  paper  presents  an  analytical  model  to  create  a  visibility  report  between  a 
range-restricted  ground  station  and  a  space  object  based  upon  the  method  of  ra¬ 
tios.  The  method  easily  lends  itself  to  implementation  on  a  personal  computer, 
processing  thousands  of  objects  spanning  the  range  of  closed  orbits  against  any  site 
latitude.  The  method  of  ratios  can  also  be  used  to  perform  parametric  studies; 
identifying  the  payoff  of  various  upgrade  strategies  for  a  tracking  station  net¬ 
work.  Figure  5  shows  an  example  of  this  application  by  varying  site  latitude  and 
sensor  range  using  mean  minutes  per  day  as  the  performance  index.  Whether 
attempting  to  quantify  the  cost  per  unit  increase  in  ground  system  performance 
or  allocating  funds  to  improve  or  maintain  a  network  of  stations,  the  merits  of 
this  tool  are  easily  recognized. 

Also,  the  method  of  ratios  can  be  used  to  identify  the  optimal  site  latitude  to 
track  a  class  of  orbits.  Figure  6  illustrates  this  process  by  holding  sensor  range 
fixed  while  varying  the  site  latitude.  Depending  on  the  satellite  viewing  priority, 
this  tool  can  help  mission  designers  position  a  ground  station  to  maximize  the 
mean  cumulative  visibility  or  reposition  mobile  tracking  terminals  to  maximize 
contact  once  a  space  vehicle  has  maneuvered.  This  method  can  also  be  used  in 
selecting  orbital  characteristics  to  maximize  contact  with  a  range-restricted  site, 
supporting  space-mission  designs  for  highly  maneuverable  vehicles. 


Mean  Visibility  (min/day) 


The  Method  ot  Ratios 


T  Range  =  1 000  km  — *—  Range  =  2000  km  — —  Range  =  4000  km 


H<j.  5.  Visibility  Report  for  Space  Object  #9  versus  Varying  Site  Latitude  and  Sensor  Range. 

Appendix  A 

Algorithm  for  the  Method  of  Ratios: 

1.  Set  L  equal  to  the  site  latitude. 

2.  Set  L'  equal  to  the  object  sublatitudc. 

3.  Set  R  equal  to  object  geocentric  range. 

4.  If  R  cos (L  -  L  )  <  1.0,  then  stop;  object  is  below  site  horizon,  0  =  0. 

5.  If  {R:  +  1  -  2Rcos{L  -  L')}  >  l> im  ,  then  stop;  object  is  beyond  site 
range,  0  =  0. 

6.  C  ompute  the  numerator  of  the  arc  cosine  argument.  If  pmax  >  {Rl  -  1). 
then 

NUM  =  \/R  -  sin(L)sin(L') 

else 

NUM  =  — — —  —  sin(L)sin(L'). 

2r\ 

1.  Compute  the  denominator,  DEN  =  cos(L)cos(L  ). 

8.  If  DEN  is  less  than  some  small  tolerance  or  if  \NUM/DEN\  >  1,  then 


Mean  Visibility  ( 
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Space  Object  #8  —  Space  Object  #9 


FIG.  6.  Visibility  Report  for  Space  Objects  and  #>)  versus  Varying  Site  Latitude  and 

4<XKI  km  Sensor  Range. 


Appendix  B 

Algorithm  for  the  Visibility  Criterion: 

1.  Get  8mAX  from  ratio  equation. 

2.  If  0MAX  is  2 it,  then  stop,  criterion  is  satisfied. 

3.  Compute  ft®. 

4.  Compute  N,uty  and  Nn<f. 

5.  Compute  and  evaluate  x,  for  an  integer  value  letting/  range  from  1  to  Nn*. 


The  Method  of  R&;io$ 
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